Proposal to stabilize and detect half-quantum vortices in strontium ruthenate thin 
films: Non-Abelian braiding statistics of vortex matter in a + ipy superconductor. 
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We propose a simple way to stabilize half-quantum vortices in superconducting strontium ruthen- 
ate, assuming the order parameter is of chiral px + ipy symmetry, as is suggested by recent ex- 
periments. The method, first given by Salomaa and Volovik in the context of HeUum-3, is very 
naturally suited for strontium ruthenate, which has a layered, quasi-two-dimensional, perovskite 
crystal structure. We propose possible experiments to detect their non abefian-braiding statistics. 
These experiments are of potential importance for topological quantum computation. 
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Introduction. The realizable prospect for topologi- 
cal quantum computation using fractional quantum Hall 
states in two-dimensional {2D) semiconductor structures 
has been discussed recently in the literature 
In particular, non-Abelian braiding statistics, manipu- 
lated by braiding particles as desired, is the robust quan- 
tum mechanical resource underlying quantum computa- 
tion in these systems. Such non-Abelian quasiparticle 
statistics arise naturally in certain exotic classes of 2D 
fractional quantum Hall states (e.g. the so-called Pfaf- 
fian quantum Hall state which may exist at the axy = § x 
plateau giving rise to the recent excitement 

about the possibility of quasiparticle manipulation and 
braiding (and eventually to topological quantum compu- 
tation) in high-mobility semiconductor heterostructures. 
Given the great interest in non-Abelian statistics and 
topological quantum computation, it is relevant to ask 
whether other physical systems could be studied exper- 
imentally to investigate non-Abelian braiding statistics. 
In this Letter, we suggest quantum vortex matter in su- 
perconducting strontium ruthenate as a prospective can- 
didate for non-Abelian statistics and topological quan- 
tum computation. The physical basis for this proposal is 
the possible px+ipy symmetry of the superconducting or- 
der parameter in strontium ruthenate. This order param- 
eter symmetry admits half-quantum vortices, which are 
analogous to the fractionally-charged non-Abelian quasi- 
particles of the Pfafhan fractional quantum Hall state. 

Superconductivity in Sr2RuOi. Spin triplet (5" — 1), 
odd parity (L = 1) superfluidity is realized in Helium- 
3 0. Recently, a body of evidence, including sensitiv- 
ity to impurity , nuclear magnetic resonance (NMR) 
Knight shift jsl, and polarized neutron scattering |3l 
have strongly indicated that spin-triplet superconduc- 
tivity is also realized in strontium ruthenate (Sr2Ru04) 
jl2l |. Furthermore, within the superconducting phase of 
this metallic oxide, time-reversal symmetry breaking has 
been confirmed by muon spin rotation 'lCr|. The exis- 



tence of a two-component order parameter is suggested 
by small angle neutron scattering exploring the flux lat- 
tice in the mixed state 11,]. Taken together, these ex- 
periments point to the conclusion |l2l Il3 | that strontium 
ruthenate, with a highest of about 1.5 K, realizes the 
simplest unitary T-breaking p-wave state allowed in its 
tetragonal crystal structure, namely px -\- ipy, which is 
analogous to the A-phase of He 3. Very recently, phase 
sensitive superconducting quantum interference device 
(SQUID) measurements 'l4| have given almost conclusive 
evidence of p-wave superconductivity 13] in Sr2Ru04. 

For spin triplet superconductors, the order parameter 
is a matrix in spin space 0, 



a/3 ■ 



(1) 



Here a, (3 are spin indices, k is the wave vector, cri.2,3 
are the Pauli matrices, and d{k) is a complex 3-vector. 
Orbital p-wave symmetry implies d^ik) — 
The tensor field d^j is the appropriate order parameter. 
For He 3 A-phase, or px -f ipy superconductor. 



= Ao{T)df_,{mj + in.j)e 



icj} 
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where Ao(r) is the temperature-dependent magnitude of 
the order parameter, m and h are mutually orthogonal 
unit vectors in orbital space, and 4> is the phase angle. 
Physically, d is the unit vector in spin space on which the 
projection of the Cooper pair spin is zero, and I — m x n 
is a preferred direction in the orbital space giving the 
direction of the Cooper pair angular momentum. 

For topological classification of vortices 16], note that 
the full symmetry group G=S0(3)x S0(3)x U(l), im- 
plying rotational symmetries in the spin and the orbital 
spaces, and the global gauge symmetry. The Px + ipy 
state has residual symmetry H=U(l)x U(l)x Z2, where 
the first U(l) is for spin rotation symmetry about the 
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vector d, the second U(l) is the combined gauge-orbital 
symmetry, and the Z2, most important for our purposes, 
signifies a discrete combined gauge-spin rotation symme- 
try of the order parameter, {d,(j)) {—d,(j) + tt). 
The first homotopy group of the degeneracy space R then 
corresponds to 7ri(R)=7ri(G/H)= 7ri((S2xS03)/Z2)=Z4, 
the cychc group with 4 elements. Here denotes the 
surface of a sphere. The vortices, then, correspond to 
only four classes, with vorticities TV — —^,0, 5, 1 — 1. 
The surprising fact is the existence of N — ±i, which 
correspond to half flux quantum per vortex. The half- 
quantum vortex is crucial for non-Abelian statistics in 
this system. 

Stabilizing Half- Quantum Vortices. The simplest rep- 
resentation of the class = ±^ is given by, 

d = xcos(0/2) +ysin(6l/2), = ±6'/2, 

/ = const. (3) 

Here, it has been assumed, without any loss of generality, 
that d lies in the x — y plane and 6 is the polar angle 
with the vortex core as the origin. The flux enclosed 
by a half- vortex is half the flux quantum, ^^o, where 
(^0 = ff • The existence of the half- vortex is possible 
because d can change sign upon circulating the core, and 
the phase angle (j) can simultaneously change by tt to keep 
the order parameter, Eq. [21 unchanged. 

However, in real systems, half-vortices are energetically 
costlier than vortices with = 1, 2 ( a vortex with N = 2 
is continuously deformable to iV = 0). This is because 
of small, but non-zero, spin-orbit or dipole energy, Eso = 
—flsoid-l)"^ , which binds I and d parallel, and so penalizes 
spatial modulation of d around the vortex core. The 
texture for a half-quantum vortex requires d and I to be 
nonparallel over a region of macroscopic size. Because of 
i?soi this region forms a domain wall terminated by the 
vortex, expending energy proportional to the volume of 
the domain wall. Another way of looking at this is to note 
that in the case of dipole- locking of d and I, the extra Z2 
symmetry of H is absent. The first homotopy group of 
R is then 7ri(S^xS03)=Z2. There are only two classes of 
vortices, the uniform class, A^ = (a vortex with N = 2 
remains continuously deformable to this class), and the 
class with vorticity N = 1. Half-quantum vortices are 
no longer possible. Hence, to stabilize the half-quantum 
vortices, the spin-orbit energy must be neutralized. 

A simple way to neutralize the spin-orbit energy, first 
proposed in the context of He 3 [ig, is very naturally 
suited to strontium ruthenate. In He 3 A-phase, one has 
to use an artificial parallel plate geometry so that the di- 
rection of the Cooper pair angular momentum, I, is con- 
strained to be along the plate normals, n. In the presence 
of a sufhciently strong magnetic field H \\ n, the Cooper 
pair spins will align themselves with H. The vector d 
is then in the plane of the parallel plates, expending the 
maximum spin-orbit energy. However, this energy cost is 



a constant for all types of vortices, d is also free to ro- 
tate in the plane to —d spending no additional spin-orbit 
energy. In this way one can avoid the large domain wall 
energy cost for the half-quantum vortices making them 
competitive with other topological excitations. 

However, this restrictive geometry is completely natu- 
ral for strontium ruthenate. It is a metallic oxide with 
highly two-dimensional (2D), layered, tetragonal crystal 
structure. In fact, the crystal structure of Sr2Ru04 is al- 
most identical to that of perovskite (La, Sr)2Cu04 high 
temperature superconductor Two dimensional, cor- 
related Ru02 layers are separated by intermediate stron- 
tium reservoir. Electronic conduction is mediated by d- 
orbitals of ruthenium ions in 2D planes. This highly pla- 
nar structure is also evident by the superconducting state 
properties such as the upper critical field Hc2, which is 
about 20 times smaller when applied parallel to the c-axis 
than when applied parallel to the ab plane. 

Sr2Ru04 thus naturally realizes a, Px + ipy supercon- 
ductor in a parallel plate geometry. It is believed that I 
is in the c-direction, and that d is pinned to / due tospin- 
orbit energy, although conclusive proof is lacking [iSj . If 
this is indeed the case, the most direct way to stabilize 
half-quantum vortices is to simply apply a magnetic field, 
strong enough to overcome the spin-orbit energy, parallel 
to the c-axis. In the case of He 3, Ref. claimed that at 
low enough temperatures T <C Tc, vortices with A^ = ^ 
are energetically slightly more favorable than those with 
A^ = 1, 2. In the absence of a definitive proof of even the 
symmetry of the order parameter being {px + ipy)-type 
and lack of knowledge of other system parameters, such 
microscopic energy comparisons for Sr2Ru04 are prema- 
ture. However, half-quantum vortices are expected to be 
at least competitive in energy and may even be more fa- 
vorable in the presence of magnetic fields simply because 
of smaller vorticity. 

The required magnetic field H can be estimated from 
the spin-orbit or dipole coupling constant flgo- Unfortu- 
nately, we are unaware of any reliable theoretical estimate 
of this quantity. In He 3 A-phase, the required H was 
estimated to be greater than ~ 25 G Sr2Ru04, a 

recent experimental estimate of the spin-orbit decoupling 
field is 250 G 0. We suggest applying a magnetic 
field of this order of magnitude in order to decouple d 
from I. Note that this can be done without destroying 
superconductivity since the upper critical field in the c- 
direction is still larger: H^2\\c ^ 750 G 0. 

Non-Abelian Statistics of Half-Quantum Vortices. 
Once stabilized, the half-quantum vortices may be de- 
tectable by exploiting their statistics under interchange 
of positions. In the following, we recapitulate the ar- 
guments leading to their non-abelian braiding statistics 
j2d I211 1231 ■ An interpretation of the half-quantum vor- 
tex can be given, when the subsystems of up and down 
spins are uncoupled, in terms of a single-quantum vor- 
tex for only one of the spins and zero vorticity for the 
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other spin component |^ |2^. This can be seen by 
explicitly writing out the components of the order pa- 
rameter matrix, Eq. ^ Ignoring, for the moment, the 
orbital dependencies of all quantities, we get, A|| — 
AQexp{i(j)){—dx +idy), Aj^^ — Aoeyi'p(i4)){dx +idy) and 
— A|| — Aq eKp(i(j3)dz. Using the form of d and (j) 
in Eq.|21 with (j) = 9/2 (positive half-quantum vortex), it 
is easy to see that, A^^ = — Aq, A^ = Aoexp(z6') and 
A-fj^ = A|-f = 0. Because the phase-angle cj) = 9 changes 
by 27r upon circulating the core, the down spin pairs see a 
full single quantum vortex, while the up spin pairs see no 
vortex at all. Because of this property, the Bogoliubov- 
de-Gennes equations separate into two pieces: the part 
for the down spin component is identical to that for an 
ordinary single-quantum vortex in a p-wave superconduc- 
tor (2^, albeit for spinless electrons, while the part for 
the up spin component is devoid of any vortex. As in 
Ref. |24| , the low-energy spectrum in the vortex core for 

the down spins is given by £"„ = nujo, where ujq = 
with tp the Fermi energy, and n an integer. For n = 0, 
there is a zero energy state in the vortex core. 

It is interesting to note that the Bogoliubov quasipar- 
ticles corresponding to the levels En are linear super- 
positions of particle and hole operators from the same 
species, namely, down spin electrons, 7I = uipl^ -f v^ii, 
where ipii is the electron annihilation field for down 
spins at the location of the i-th vortex. They satisfy 
7|(i?„) = 7j(— £"„). The zero energy state, then, is a Ma- 
jorana fermion state, satisfying 7/(0) = 7j(0). Each half- 
quantum vortex, i, has a zero energy Majorana fermion 7; 
in its core. In the case of 2m such spatially separated vor- 
tices, the 2m Majorana operators can be pairwise com- 
bined to produce m complex fermion operators, each of 
which can be occupied or unoccupied in the ground state. 
The ground state, thus, is 2™ fold degenerate. Further, 
using the property that 7i's carry odd charge with re- 
spect to the gauge field of a single-quantum vortex, that 
is, 7,; — > — 7i for a phase change of 2ti, it follows that 
upon interchange of two neighboring half-quantum 
vortices 7^, 7^+1, 

7j -> 7^+1, 
7,+i -7j 

1] ^ I3 j T^hj ¥"1 + ^ (4) 

Observing non-Ahelian Braiding Statistics. These 
topological properties are identical to those of the Pfaf- 
fian quantum Hall state, which is a candidate descrip- 
tion of the observed fractional quantum Hall plateau with 
'^^v ~ 1 1[ '^"^ two-dimensional electron gases 
This state is distinguished by the non-abelian braiding 
statistics of its quasiparticles, which are precisely 
This state has recently been the subject of much inter- 
est as a platform for fault-tolerant quantum computation 
. Recent papers B, Q , elaborating on an earlier pro- 
posed experiment j2y| explain how non-Abelian statistics 



could be directly detected through quantum interference 
measurements. In this paper, we suggest that the same 
underlying physics could arise in Sr2Ru04. Following ref. 
Q, we note that the two possible states of a pair of vor- 
tices in a 2D px + ipy superconductor (corresponding to 
the presence or absence of a neutral fermion) can be used 
as the two states of a qubit. When the two vortices are 
far apart, no local measurement can distinguish the two 
states of the qubit, so they form a decoherence-free sub- 
space. However, there are important differences between 
Sr2Ru04 and 2DEGs which lead to practical differences, 
which we discuss below. 

Vortices in Sr2Ru04 are neutral, unlike quasiparticles 
in a quantum Hall liquid. Therefore, a direct quantum 
interference measurement of quasiparticle transport as 
proposed in ^ would not apply to Sr2Ru04. However, 
the Josephson relation tells us that the motion of a vortex 
across a line connecting points 1 and 2 causes a phase 
slip, thereby inducing a voltage drop between 1 and 2 
according to 

K = -i™ (5) 

e " 

When vortices become depinned, either as a result of non- 
equilibrium effects or thermal fluctuations, there is dissi- 
pation. 

A second important difference is that, even in a thin 
film, there are many 2D layers in a Sr2Ru04 crystal. 
Hence, a vortex which runs through N layers actually 
has N Majorana zero-modes associated with it, if we ne- 
glect inter layer hopping. Consequently, a pair of vortices 
has N qubits associated with it. When two vortices are 
braided, their topological interactions, as summarized in 
21 occur independently in each layer. In other words, the 
Majorana modes in layer a only affect other Majorana 
modes in layer a, not those in layer /3. Interlayer hopping 
causes these fermionic zero-modes to hybridize and form 
a band along the c-axis. While these states are split in 
energy, their topological protection is unaffected. 

Another important difference is that vortices are po- 
tentially heavy objects, so it is less clear that one can ob- 
serve the coherent quantum motion of vortices. However, 
vortices are expected to become light as the Mott transi- 
tion is approached. While the normal state of Sr2Ru04 
appears to be a Fermi liquid at low temperatures^Jrans- 
port at temperatures above 30 K is anomalous 13]. If 
this is an indication of proximity to the Mott transition, 
then it could imply that vortices are light. In general, 
we expect vortex motion to be coherent in the absence 
of appreciable quasiparticle excitations. Having very low 
temperature would help prevent decoherence by quasi- 
particles and by phonons. 

Keeping all of these caveats in mind, we turn to exper- 
iments which could observe the topological properties of 
vortices in Sr2Ru04 thereby, in principle, taking the first 
steps towards topologically protected qubits. Our goal 
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is to observe the two topologically-degenerate states of a 
vortex. 

Recently, Stone and Chung emphasized that when 
two vortices are brought together, the difference between 
the two states of a qubit is now manifest: there either 
is or isn't a neutral fermion in the core of the combined 
vortex. There is an energy cost in the latter case. Since 
we expect neither an excess of charge nor spin in the vor- 
tex core, we can look for subtle differences in the charge 
distribution or try to observe the energy difference. The 
idea is the following. While the total charge is the same 
whether or not the neutral fermion is present, the lo- 
cal charge distribution will be different in the two cases. 
Thus, if one had spatial resolution smaller than the size 
of a vortex core ~ 660 A, one could try to look with 
an STM for the particular charge density profiles associ- 
ated with the presence or absence of a neutral fermion. 
In principle, gravitational effects can be used to detect 
this energy difference, but the required sensitivity may 
be very difficult to achieve. 

An alternate strategy is to observe the braiding statis- 
tics of vortices through their effect on quantum interfer- 
ence. Consider the experimental setup in Figure We 
assume that a thin film has a hole bored in it so that it 
has an annular topology. Flux $ is threaded through the 
hole, so that there are iV^ vortices in the hole. A mag- 
netic field H > Hci is applied to the rest of the film so 
that vortices penetrate the rest of the film as well. A cur- 
rent J is driven in the y-direction. Vortex motion in the 
x-direction generates a voltage drop in the y-direction, as 
in ©. If the motion of a vortex is coherent, then the two 
trajectories in the figure will contribute to the resistivity 
according to |2^ : 

Ryy (X \ti\^ + \t2\^ + 2Rc{tit; (6) 

where ti, ^2 are the amplitudes for the two trajectories, 
B is the braiding operator for a vortex to go around the 
vortices in the hole, and \^) is the combined state of 
the moving vortex and the vortices in the hole. When 
the number of vortices in the hole is even, \B\ ^P) is 
a phase, e*^. However, when the number of vortices in 
the hole is odd, \B\ *) 0, If we can vary ti 
and t2, then we can distinguish these two cases. It is 
not clear how to do this, but there are several possibili- 
ties. (1) If we gate the Sr2Ru04 thin film along one of 
the possible trajectories and apply a gate voltage, this 
will suppress the electron density along this trajectory. 
This, in turn, will change the amplitude for this trajec- 
tory. (2) If pressure is applied to the part of the film 
where one of the trajectories is located, the lattice con- 
stant will be different there, presumably translating into 
a different amplitude for this trajectory. (3) Applying 
a temperature gradient will cause the superfiuid density 
to be different along the two different trajectories which, 
again, would lead to a variation in their respective am- 
plitudes. However, there is a caveat to this final option: 



raising the temperature may be dangerous since we want 
vortex motion to be coherent. One could also try to ob- 
serve interference in the Nernst effect in the geometry of 
Figure n but this is difficult experimentally. 
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FIG. 1: In the experimental geometry depicted above, a thin 
film of Sr2Ru04 has a central hole with a flux $ through 
it. When an electrical current flows in the y-direction, de- 
pinned vortices flow in the a;-direction. This generates a volt- 
age drop in the x direction. If the motion of a vortex is 
coherent, then there will be interference between trajectories 
which go around either side of the hole. The nature of this 
interference will depending strongly on whether the number 
of half-quantum vortices in the hole is odd or even. 

Conclusion. In summary, we propose an experimen- 
tally feasible technique for stabilizing half-quantum vor- 
tices in Sr2Ru04, assuming the system to be a p^, -I- ipy 
superconductor. If stabilized, such half-quantum vortices 
should exhibit non-Abclian braiding statistics and, po- 
tentially, lead to topological quantum computation. 
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